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Abstract 
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We study magneto-optical properties of monolayer graphene by 
£T) ■ means of quantum field theory methods in the framework of the Dirac 

model. We reveal a good agreement between the Dirac model and a re- 
—l. | cent experiment on giant Faraday rotation in cyclotron resonance [19]. 

We also predict other regimes when the effects are well pronounced. 
(^ , The general dependence of the Faraday rotation and absorption on 

CN ' various parameters of samples is revealed both for suspended and epi- 

taxial graphene. 

>: 
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u ■ 1 Introduction 

a ■ 

Graphene, which is a one-atom thick layer of carbon atoms, does not need 
any lengthy introduction nowadays. Ever since its experimental discov- 
ery, graphene is in the spotlight of applied research in condensed matter 
physics [1, 2, 3, 4]. The dynamics of quasiparticles in graphene is governed 
by the Dirac model [5], i.e. they satisfy the Dirac equation where the speed 



of light c is replaced by the so-called Fermi velocity vp — c/300. Due to 
this quasi-relativistic nature of electronic excitations, much of the graphene 
physics is in fact described by Quantum Field Theory (QFT), see [6, 7] for 
recent reviews. 

One the main obstacles to further industrialization of graphene is the 
gapless spectrum of quasi-excitations, and thus most of the recent research 
was devoted to investigation of the electron transport in graphene sheets 
and possible ways of generating a controllable mass (energy) gap, see discus- 
sion in [4]. On the other hand, the optical properties of monolayer graphene 
from the beginning proved themselves extraordinary The very first measure- 
ment [8] of optical absorption in graphene revealed a surprisingly huge effect 
- monoatomic layers absorbed about 2.3% of the incident light, exactly as 
predicted by the Dirac model, see [9, 10, 11, 12]. 

In this paper we investigate the optical properties of monolayer graphene 
samples in external magnetic field. We focus here on the Faraday effect. It 
may be thought of as an analog of the Hall effect considered at "non-zero 
frequencies" , and this similarity was used in Ref. [13] to conjecture that the 
former should be common for Hall systems. In Ref. [14] the Faraday rotation 
was related to possible non-compensation of parity-odd parts of the polariza- 
tion tensor between various generations of fermions. Other related works on 
magneto-optical properties of monolayer graphene include theoretical studies 
of the Kiev group, e.g. [30, 32], and numerical simulations in [15]. 

The quantity and quality of magneto-optical experimental data is grow- 
ing. In the case of multilayer graphene samples the first measurements of 
transmission spectra were reported in [16], and the work was later continued 
in a series of papers, e.g. [17, 18]. More recently, both the Faraday rotation 
and transmission in mono and multilayer epitaxial graphene were measured 
in [19]. The Faraday rotation appeared to be unexpectedly large and was 
dubbed the "giant Faraday effect". The experimental setup of [19] permitted 
a detailed investigation of the so called cyclotron resonance regime. 

To explore these effects we use a quantum field theory approach to the 
physics of graphene. Quasi-relativistic dispersion relation of the electronic 
excitations in graphene makes QFT a more adequate language to describe its 
properties as opposed to non-relativistic quantum mechanics. Application of 
the field-theoretical methods allow one to investigate such purely relativistic 
effects as the Klein paradox, or appearance of gauge fields in graphene. It 
also permits to investigate other ideas of QFT and quantum gravity in low 
energies bench-top experiments - - a situation hardly imaginable even ten 



years ago. 

Following this approach, we start with the Dirac-Maxwell action, where 
the spinors are confined to the graphene surface, while the electromagnetic 
field lives in the ambient 3 + 1 dimensional space. Electromagnetic field is 
not quantized, and is given by a sum of a constant magnetic field and a 
fluctuating one (photons). The magnetic field is taken into account exactly, 
while in photons we keep the second order terms only. Therefore, on the 
quantum field theory side we calculate a single one loop diagram (though 
a quite complicated one, as the propagators contain the magnetic field, the 
chemical potential, the mass, and a phenomenological parameter which de- 
scribes impurities). This calculation is similar to that of [30, 32], but is done 
differently to facilitate a comparison to experimental data. 

With some assumptions on the parameters (frequency, magnetic field, 
impurities, etc) one can derive expression for conductivity of graphene and 
similar materials [38, 39] that are simpler than the ones given below. We like 
to stress, however, that the conductivity that we calculate and use is an exact 
result of quantized Dirac model at the first order of the fine-structure con- 
stant a but without any restrictions on the other parameters. The advantages 
of having such expressions (perhaps, complicated) are obvious. Among them 
is the possibility of a cleaner judgement of validity of this or that approxi- 
mation. The price to pay is the use of the QFT machinery and terminology. 

While the experiment of [19] can be fitted by the Drude model, the use 
of complete Dirac model makes it possible to relate the optical properties of 
graphene to its microscopic parameters. It also permits to investigate and 
predict other regimes (besides the cyclotron resonance one studied in [19]), 
where both the Faraday rotation and absorption are well pronounced. In 
particular, we reveal quantization of the Faraday angle for clean samples, 
and predict pole-like peaks for relatively high frequencies. 

Similar results were recently obtained in [20] using an approximate method 
of 'equation of motion' for calculation of the conductivities of graphene. 

We also like to mention the theoretical studies [40, 41] of the Faraday 
effect and related phenomena in topological insulators, which are closely 
related to graphene [4]. 

The paper is organized as following. In the next section we shortly review 
the basics of the Dirac model in graphene. In Section 3 we formulate its 
connection to optical properties of both epitaxial and suspended samples. In 
Section 4 we present the calculation of the polarization operator in magnetic 
field both for ideally clean and realistic samples, and finally in Section 5 we 



investigate the Faraday rotation and transmission of light performing the fit 
of existing experimental data, and predicting new regimes. 

2 The Dirac Model 

The Dirac model of the quasi-particles in graphene is based on the tight- 
binding model 1 . An interested reader can find its detailed description in a 
number of reviews, e.g. [5, 31], while we briefly formulate the model itself. 

Consider an (infinite) graphene sheet occupying the plane x 3 = in a 
3 + 1 dimensional ambient space. The action for an electromagnetic field 
interacting with quasiparticles propagating in the graphene surface reads 

S=-\j d 4 x F^ + J d 3 xi>$il> (1) 

with F^ = d^A v - d v A^ 

Jj) = Jj){A)= if(dj + ieAj) - A . (2) 

Greek letters denote the coordinates of the ambient space, /j,,is — 0,1,2,3. 

Latin letters j,k = 0, 1, 2 correspond to coordinates on the graphene plane. 

Since there are N = 4 species of fermions in graphene, the gamma matrices 

are in fact 8x8, being a direct sum of four 2x2 representations (with two 

copies of each of the two inequivalent ones). The matrices 7 J are rescaled, 

7° = 7°, 7 1 ' 2 = vp'-f 1,2 , as compared to ordinary gamma-matrices 7q = 

— (7 1 ' 2 ) 2 = 1. v f — 1/(300) is the Fermi velocity. We use the units where 

H — c — 1. A is the mass gap. The Maxwell action is normalized in such a 

way that 

Atx 
e 2 = 4vra = — . (3) 

By averaging over the fermions in (1) one arrives at an effective action 
for the electromagnetic field, so that the full action becomes 

S A = - 1 -Jd'xF^ + S cS (A) (4) 



1 Magneto-optical properties of graphene can be analyzed [42, 43] starting directly with 
the tight-binding model. This is technically more complicated than the analysis based on 
the Dirac model. Both approaches coincide for the energies below ~ 4.4eV [43]. 



where 

S eS (A) = -ihidet(lp(A)). (5) 

Since the action (1) is quadratic in ip the expression (4) is exact. 

We shall be interested in the scattering of light by a graphene sheet placed 
in a strong external electromagnetic field. Therefore, it is convenient to split 
the electromagnetic potential into a background and fluctuating parts, 

A, -► A, + A*. (6) 

In the case of the Faraday effect the background part A b describes a con- 
stant magnetic field B perpendicular to the graphene plate, while A^ is the 
photons. The background part has to be taken into account exactly. On the 
other hand, for our purposes it is enough to expand S e s up to the quadratic 
order in the radiation field A^, which is given by the one-loop diagram 

S cS (A) = A-^y^A 

= U^ Aj{ - p)uJk{p;Ab)Mp) ' (7) 

where the polarization operator ELy* 
W k (p) 
= *e 2 / ^^tr [S(q , q)f 'S(q -p ,q- p)^ k ] (8) 



depends on A b through the fermion propagator S = {Ip{A )) -1 . Here and 
below with bold Latin letters we denote spatial in-plane vectors, p = (pi,p2)- 
The propagation of electromagnetic waves in the whole space can be de- 
scribed now by modified Maxwell equations with a delta-function interaction 
corresponding to (7) 

dpF*" + S{x 3 )U up A p = , (9) 

where we set IT 3At = IP 3 = 0, p = 0,1,2,3. These equations describe free 
propagation of light away of the graphene subject to the matching conditions 

4 1 — A I 

Slfl\ x 3 =+0 — Si[j,\ x 3 = _ /^px 

(dsA^+o - (^A A1 ) S 3=_o = n„ »A u \ x3=0 {iU) 



on its surface. The same delta- function term in the Maxwell equations can 
also be interpreted as an in-plane current. Then the ac conductivity, that 
is a proportionality coefficient between the current j and electric field E, 
ja = o'abEb, is expressed as 

n ab (w) 

<?ab = : 

where a, b = 1, 2. This relation, where the polarization operator is considered 
independent of the spatial in-plane momenta, i.e. describing the normal 
incidence case, can be inverted 

W k (u,p = 0) = 

^ tu (5 ab a xx {u) + e ab a xy {u)) J' { } 

Here S ab and e ab are two-dimensional Kronecker and Levi-Civita symbols, 
respectively. By calculating explicitly the corresponding components of (8) 
we shall find both diagonal and Hall ac conductivities, see section 4. The 
temporal component, ir°(u), will be of no importance in our further consid- 
erations. 



3 Optical properties of graphene 

Let us consider a beam of light linearly polarized along the x = x l axis 
passing normally through a graphene layer on a substrate of a finite thickness 
d. In a magnetic field the Hall conductivity of graphene a xy is non-zero. 
Therefore both x- and y-polarizations will be present in the transmitted 
light. The intensity of light measured after passing graphene, substrate and 
a polarizer rotated by angle in respect to the x-axis is given by [21] 

If | 2 _i_ If | 2 |f | 2 _ |f | 2 

T \^xx\ T^ \ L xy\ \ v xx\ \ v xy\ ,-. , 

h = g 2 ^ 

-Re(t_t* )sin20, (12) 



where t XXjXy are the transmission coefficient of x, y-polarizations, respectively. 
The angle 9 where (12) attains its maximum as a function of gives the 
angle of polarization rotation 

^-larctan^H. (13) 

£• \ h xx\ \ h xy\ 



The total transmission is just a sum of the intensities in two polarizations, 

By solving (10) with II defined by (11) along with a standard matching 
condition at the dielectric-vacuum interface at z = x 3 = d, we obtain 

t xx = -±n s e^-^-^- (15) 

t xy = ln s e^-^-^- (16) 



where 



A x = (n s - l)(n s - 1 - a xx )e 2i(kvn ° - {n s + l)(n s + 1 + a xx ) 

A y = a xy (n s + l + (n s - l) e 2id " n °) , (17) 



and n s is the refractive index of (non-absorbing) substrate. 

Both polarization rotation angle and transmission oscillate rapidly with 
the frequency u (or with the substrate width d). This effect is called the 
Fabry-Perot oscillations. However, in many cases the oscillations are smeared 
out due to low resolution of the measurements, or other sources of incoher- 
ence. The averaged intensity can be obtained by integrating (12) in d over 
the period of oscillations, d p = - J! -. 

Thus, we have 

1 8 Re{a xx (T* xy ) + 4{n 2 s + 3)Rea xy 

U — ctrctjctll 



2 ' \2a xx + nl + ^-A\a xy \^-{nl-lf 

~ |a+| 2 -|M 2 + |a_| 2 -|M 2 ( ' 

where a± = (n s + l)(n s + l + a xx ±ia xy ), b± = (n s -l)(n s -l- (a xx ±icr xy )), 
T was first obtained in [17]. Without graphene, the transmission of the bare 
substrate is 

To = ^. (19) 

1 + nj 

Later on we shall use the normalized transmission, 



T = f/T , (20) 



which characterizes better the properties of graphene itself. 

To obtain the first order expressions in a one keeps only the linear terms 
in the conductivity 

T ~ 1 - 2f s -Re<j xx , ~ -f s Rea xy (21) 

where the coefficient f s reads 

f " ^^- (22) 

To reproduce the case of suspended graphene one puts n s — 1 (i.e. / s = 
1/2) to obtain 2 

/t-, , | ^xo: ~r ^ I I &xy \ 



a 2 j_2 _|_ Afl. + 4|2' 
I Z ~r Ojj T" 1(J X y 



e =-^ <::_ : ( 23 ) 

while in the linear order in the conductivity the expressions are 

T ~ 1 - Rea xx , 9 ~ --Rea xy . (24) 

4 Calculation of the polarization operator 

The polarization operator IL^ has been considered in a number of papers for 
systems characterized by various sets of parameters. In the zero-temperature 
case it was first calculated in [23]. At finite temperature calculations were 
done in [25], while some other cases were considered in, e.g., [26, 33, 35]. 
In the presence of an external magnetic field extensive calculations of the 
polarization operator were done by the Kiev group and collaborators [28, 29, 
30, 32, 33, 34]. Here we shall only sketch the derivation. 

In a constant magnetic field B perpendicular to the surface of graphene 
the fermion Green function has the form [24, 29] 

cf, , n \ v-p 2 /\eB\ VV i\ n t ^»( a; + / i 'P) / c\ 



2 The reflection and transmission coefficients for suspended graphene can be read off 
from [14] after the identification a xx = iaty /lu, <j xy = a<p. 



where 

2p' 



'2p 2 



P_L 



n I Dl 



— P+L n -i 



l ''' PTL --i(^j) • ( 26 ) 



Je-B| 
Here P± = (1 ± i r y lr y 2 €B)/2, e# = sign (5), the Landau levels are M, 



2nvp\eB\ + A 2 , and L^ are the associated Laguerre polynomials, L n = 

l4°\ i i = 0. /i denotes the Fermi energy shift (or the chemical potential 
in the QFT terminology). 

Next, we need to introduce a phenomeno logical parameter V which de- 
scribes the presence of impurities. This can be done by means of the substi- 
tution 

u> — > u> + iT sgncu, T > 0. 

everywhere in (25). In the limit r — > one recovers the usual Feynman 
propagator. For r ^ the propagator S is not an analytic function of uj 
due to the presence of sgnw. Most generally, V can dependent on frequency, 
magnetic field, etc. We will restrict ourselves to a constant T, which is 
sufficient in a not too strong magnetic field [44], but see also [45]. 

Describing the disorder in such a simplified manner, we assume that the 
long-range impurities present in considered samples of graphene are suffi- 
ciently weak (if any). Otherwise the states near the Dirac points get local- 
ized, and deviations from the Dirac dispersion should be taken into account. 
The behavior of graphene in the presence of strong long-range impurities is 
studied in detail in [46]. 

We are ready to substitute (25) in Eq. (8). Since we shall consider normal 
incidence wave only, p = in (8), the loop integral over spatial momenta q 
can be easily performed with the help of the orthogonality condition for the 
Laguerre polynomials 



x a e- x L ( n a \x)L^\x)dx = (n+ , a \ 



(27) 
o '" "" n\ 

Thus, we arrive at the following expression for the Hall conductivity 

e B aNLl ™ 
<7 xy {po) = 2^ G n (28) 

PO n=0 

9 



where L\ = 2v 2 F \eB\, G n = J*f^ du (/„,„+i - f n +i,n) and f nm is given by 

fnm = — [(w + iT sgna;)(p - w - ir sgn(u; - p )) + A 2 ] 

x((u + iT sgnou) 2 - M 2 )' 1 

x((w -po + ^ sgn(u- Po )) 2 - M 2 m )- 1 (29) 

Here u = u> + \i. For the diagonal conductivity, we have 

, , iaNL 2 h ~ TT 

a xx ( Po ) = b - £ H n (30) 

Po n=0 

where H n = J^ du (/ n ,„+i + / n +i,n)- 

Note, that since for describing graphene the 7-matrices are taken in a re- 
ducible representation consisting of equal number of inequivalent irreducible 
ones related to each other though the parity transformation, we have, in 
particular, tr (7°7 1 7 2 ) = 0. 

4.1 QHE in clean graphene 

Calculations that we present here are rather standard, but differ in details 
from other sources. As a consistency check we reproduce in this subsection 
the Hall effect in graphene. 

For a clean graphene, T — > +0, the calculation of the frequency integral 
in (28) can be performed with the help of the Cauchy theorem yielding 

9 Pi ~ M 2 - M 2 , , + 2A 2 , N 

°* M - ~ wNLl (ri + Vr-^T-4riwu (31) 

where n is defined in such a way that M no < /i, while M no+1 > //, i.e. 

/i 2 -A 2 



n 



2v 2 F \eB\ 



(32) 



here |_^J denotes the integer part of x, defined so that |_^J = for x < 0. As 
expected, (31) coincides with the previous results and, for instance, with the 
T ->■ limit of Eq. (12) in Ref. [30]. 
In the dc limit, po — * 0, one has 

0% = lim a xy = ae B N (1 + 2n ) . (33) 

10 



This again coincides with the results of other calculations of relativistic Hall 
conductivity [29, 49, 50, 52, 53]. The law (33), which is called the anomalous 
quantum Hall effect, or the unconventional integer quantum Hall effect, was 
checked on experiments [47, 48]. This is one of the most spectacular confir- 
mations of the Dirac model model of quasiparticles in graphene. Note, that 
this result holds for any value of the mass gap. 

In the opposite limit, p — >■ oo, (i.e., for the visible light), we get 

a xy ~ -^i. (34) 

4.2 Impact of impurities 

In realistic samples impurities are always present, so that one should keep 
T positive. In this case, the fermion Green function, S, is not an analytical 
function of u, and we cannot apply the residue theorem for a;-integration 
in (28) and (30). Still, either this integration or n-sum can be resolved 
explicitly there, and it is a matter of convenience which operation to perform 
first. For analysis of the conductivity dependence on the frequency we find 
it more suitable to integrate first over u. 

After rather elementary but cumbersome algebra the integration in the 
off-diagonal part gives 

1 gi (eM n ,e'M n+1 ) 



G = — V 
n 87r e ,fc± \P° - ( eM « + e ' M "+i) 

g 2 (eM n ,e'M n+l ) \ 

+ Po + 2 l T-(eM n + e>M n+1 ) ~ ^ "* ^ ' (35) 




A 



2 



M n+l M n ) 
(go + iT + n - M n+1 )(iT + n + M n+1 ; 
(po + *r + n - M n )(iT + n + M n ) 
A 2 \ 



(36) 



M n+1 M n J 
{po + tT + fi- M n ){tT + (j, - M n+1 ) 
(p + iT + n- M n+1 )(i£ + n - M n ) " 

The r — > limit taken in the above expression restores (31). 

11 



(37) 



In the diagonal part of conductivity, we obtain a similar result 

fn(eM n ,e!M n+1 ) 



H, 



n 




K f t^ + \Po-(eM n + e>M n+1 ) 



h 2 (eM n ,e'M n+l ) \ 

Po + 2«r - (eM n + e'M n+ i) 




A 



2 



M n+1 M n ; 
(ir + /i + M ra )(ir + /i + M ra+1 ) 

(po + iT + fi- M n )(p + iT + ft- M n+ i) 
A 2 \ 



(39) 



M n+1 M n ) 
{po + iT + fi- M n ){p + iT + fi- M n+1 ) 
(jT + ix - M n ){iV + ix - M n+1 ) 

It is clear from (35) and (38), that the off-diagonal part of the conductivity 
is odd in chemical potential, while the diagonal one is even. 

4.3 Renormalization 

The polarization operator (8) is power-counting divergent in ultraviolet. It 
indeed shows up in the diagonal part (30) since 

H n ~ — + 0(n- 3/2 ). (40) 

«» 2M n K ' K J 

and the infinite sum in n does not converge. As we can see, the divergency 
is in the imaginary part of a xx , while in the a 1 order we have 1 — T ~ Kea xx 
(21). Thus, if one is interested in the a 1 level only the renormalization 
procedure is not required. Generically it can be handled, for example, via 
the Pauli-Villars substraction at infinite mass. 

In doing so, one considers a difference between two polarization operators 
(8) taken at different masses A, A, where the latter one shall be taken to 
infinity after the loop momenta calculation. In (30) it results in substitution 
of H n (A) by H n (A) — H n (A). To calculate the asymptotics of such o xx at 
A — > oo one adds and subtracts (40) from under the n-summation and then 
applies the Abel-Plana formula to express the sums via integrals, performs 



12 



the expansion in large A, and then takes the limit. The renormalized diagonal 
conductivity obtained in such a way reads 



o 



R 



odNLl / ~ / 1 \ A 1 2r 



E \ H n + T^T- )+ T2 - 



Po \£>V~" ' 2MJ • LI 4A irLl 
i r°° dx ( 1 1 



2 Jo e 27rx -1 \M(ix) M(-ix) 



(41) 



where M(z) = + JzL 2 + A 2 . In Sect. 5.1 we will use this expression both 
for finite mass, and in the A — y limit. The later one is consistent with 
the use of Pauli-Vilars regularization since in 2 + 1 dimensions there are no 
logarithmic divergencies unlike the 3 + 1 dimensional case. 

For vanishing /i, A and B the one-loop diagram in (7) is finite, see [35]. In 
our case this diagram is divergent since we have considerably less symmetries. 
We would like to mention, that one has to be careful with power-counting 
divergent diagrams that are finite for symmetry reasons. Sometimes, such 
diagrams nevertheless require a finite renormalization. For example, this 
effect gives rise to to the parity anomaly [36, 37] in 2 + 1 dimensions. 



5 Phenomenology 

The above formulas can be used now to investigate optical properties of 
realistic graphene samples, including those on a substrate. In particular, in 
the next subsection we show that the Dirac model nicely fits the results of 
the experimental measurements of magneto-optical properties of epitaxial 
graphene [19]. 

Although the polarization rotation and absorption are very strong for a 
one-atom thick material, these effects are still sufficiently small to allow for 
a qualitative description by the a 1 -order of the perturbation theory. It is 
sufficient to use (21) and (28,30) in order to estimate the order of magnitude 
of the effect for transmission and rotation angle 

_/\f r2 oo 

PO n=0 

9 = -f s e -^^Re±G n . (42) 

Po „=o 
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Still, for the correct numerical comparison with the results of the experiment 
one should use the exact formulas (18) as discussed in the next subsection. 

To distinguish between different physical regimes, we shall start by inves- 
tigating the analytical behavior of both T and 9. Contrary to the naive 
expectations, neither G n nor H n as functions of frequency po has actu- 
ally any poles in the complex plane. Indeed, if you consider, for example, 
Qzipoi —M n , M n+ i) (36) in the vicinity of p = M n+1 — M n — 2iT you find out 
that it vanishes there and thus cancels out the corresponding denominator. 
Similarly it behaves at all other points eM n + e'M n+ i, etc. 

On the other hand, in the clean graphene the Hall conductivity (31) (and 
thus the polarization rotation angle) has poles at 

p = ±(M„ 0+1 ±M no ). (43) 

One would naturally expect that the behavior of conductivity in realistic 
graphene samples with sufficiently small F-s is somewhat similar to the clean 
case. Thus, in what follows we shall concentrate on these points and inves- 
tigate the optical properties in their vicinity. 

Two main physical regimes can be distinguished already by consider- 
ing (43). For a large chemical potential, /i — ¥ oo, the system enters into 
the so called cyclotron resonance regime, which is characteristic to epitaxial 
graphene in magnetic field. The other limit, \x — > 0, represents suspended 
graphene samples without gate voltage applied. The (positive) frequencies 
(43) in these cases are 

L 2 

Po 2 ^' Po ;r 2 ^' Po -n M i- ( 44 ) 

In what follows we consider these two regimes in more details to correctly 
estimate the magnitude of the effects, and specify the position of their max- 
ima. 



5.1 Epitaxial graphene, big chemical potential 

The epitaxial graphene is characterized by a considerable Fermi energy shift 
due to the interaction with the atoms of the substrate. In known experimen- 
tal devices it is of the order of tenths of eV, /j, = ep ~ 0.1 -r 0.5eV [17, 19, 27]. 
If a gate voltage is applied the chemical potential can be increased further 
more. The other parameters of the system such as the scattering rate or 

14 




Figure 1: Experimental data (dotted lines) and the best fit curves (solid 
lines, /i = 289.2meV, T = 4.4meV) for the Faraday rotation and relative 
transmission (insert) at B = 7T as functions of po[ me V]. 

the distance between corresponding neighboring Landau levels (at moderate 
magnetic fields of the several Tesla) are of the order of 1 -j- lOmeV, or smaller. 
Thus the chemical potential gives the highest scale, and the system is in the 
cyclotron resonance regime. 

The experimental measurements of the Faraday rotation and intensity 
transmission in this regime were recently reported in [19]. In their set-up the 
authors used mono and multilayered graphene samples on a SiC substrate 
subject to a magnetic field varying from 0.5T to 7T. It is the aim of the 
current subsection to show that the Dirac model fits nicely the experimental 
curves for monolayer samples. In what follows we focus on the Faraday 
rotation angle, but similar considerations can be given to the transmission 
as well. 

Among the two points given in (44) for large chemical potential, we shall 
consider the first one only, since the second one po — 2/x corresponds to the 
frequencies that are too high. Qualitatively, the behavior of the Faraday 
rotation angle at the second point is similar to the one depicted on Fig. 3 
below. 

Before performing the numerical fitting we note that the main contribu- 
tion to the infinite sum in (42) comes from a small vicinity of n = no([i), 
(32). It follows from analytical investigation of ReG n , which also shows that 
for big enough /i, roughly for // > ^-p, each contribution to the sum has a 
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maximum as a function of frequency at po = 

d(G n / Po ) 



Ro- 



dp 



= 0. 

Po=0 



The maximum value, ReG n (po = 0), as a function of n for fixed chemical 
potential, is sharply peaked around n = n , the localization being governed 
by the difference of two slightly dislocated arctangents. 

Thus, for comparing with observations we can replace the infinite sum- 
mation in (28,30) by a finite sum. For the given experimental data of [19], it 
proved to be sufficient to choose the summation range for n from to n + 30. 
A rough estimate of the amplitude of the effect can be given by the main 
contribution with n = n , 

_ € B aNL 2 b 

"max — Js -ti-eLTjjQ 

PO p =Q 

0.75f s aL 2 b 3f s afi 2 



(M no+1 -M„ () ) 2 L 2 ' 

which is in a reasonable agreement with direct numerical calculations (as well 
as with experimental results, provided we put f s ~ 0.31 according to [19]). 

Despite that the approximate formulas (42) give a correct order of mag- 
nitude estimate of the observed effect, to perform precise numerical fitting 
of the experimental data we should use the exact formulas (18). The results 
of the fitting which was performed simultaneously for polarization rotation 
and relative transmission (as defined by (20)), are summarized below: 

• The experimental curves for rotation and transmission are nicely fitted 
by the Dirac model at fields greater then 2T. The fit is exemplified on 
Fig. 1 where both experimental (dotted lines) and theoretical curves 
(solid ones) are presented. 

• The best fit value of the chemical potential fi (i.e. the parameter which 
enters the Dirac model) is found for such fields to be approximately 
280 -7- 290meV, somewhat lower then the Fermi energy shift measured 
at zero field, €f = 340meV. It does not show any specific dependence 
on the magnetic field. 

• The fit turns significantly worse for magnetic field of 2T and lower. 
However, it can be somewhat improved by greatly lowering the chemical 
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Figure 2: The rotation and transmission (insert) theoretical predictions (solid 
lines) plotted against the frequency po[ me V] at B= 2T for fi = 120meV, 
T = 3.9meV, A = HOmeV. Corresponding experimental curves are given in 
blue dotted lines. 

potential (down to about 100 -j- 120meV) and introducing a non-zero 
mass gap A, see Fig. 2. 

We do not perform a statistical analysis of the data and do not calculate 
the accuracy of determination of the parameters //, T and A. However, some 
more qualitative remarks are in order. 

The experiment was conducted at temperature of 5K and therefore the 
zero-temperature conductivity considered above is sufficient to describe the 
observed data. The dependence of best fit value of T on the magnetic filed 
was found to be about ±10%, and we were not able to resolve it further more 
at the given accuracy of the experimental data. Our calculations showed al- 
most no dependence on the mass gap A, provided it is much smaller the 
chemical potential fi. On the other hand, a mass gap close to the value of 
chemical potential does change the curves significantly, as on Fig. 2. How- 
ever, one should not consider such fitting procedure as a strong indication 
of a mass-generation in given graphene samples. It rather calls for more de- 
tailed investigation of the phenomena at low fields. The latest research, [22] , 
proves that the spectral shape of the Faraday rotation and the transmission 
can be affected by magneto-plasmonic effects which are of greater impor- 
tance for lower fields. These effects occur due to nanoscale inomogeneities in 
epitaxial graphene which, in principle, can be taken account by introducing 
a background field in the Dirac action. This is an interesting problem, that 
we are going to address in the future. 
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In the original paper [19], the experiment was fitted using a linear order 
approximation in conductivity, the later described by the Drude model. It 
was found there that the cyclotron frequency dependence on magnetic field 
does not conform the theoretical prediction (44), u c ~ Li/2/x, first obtained 
in [32]. 

Apart from the deviation of the cyclotron frequency from its' theoretical 
value [19], the Drude model was found to be nicely describing the experi- 
ment. Still the Drude-like behavior of the Dirac model cannot be uniformly 
extended for higher frequencies since the later does not predict any poles of 
conductivity in the complex plane showing a much more complicated behav- 
ior. 

Finally, we would like to mention here that calculations performed in the 
framework of the equation of motion method [20] are also in a good agreement 
with the experimental results of [19]. 

5.2 Suspended graphene: moderate and small chemi- 
cal potential 

We continue by investigating the Faraday rotation as a function of frequency 
for regimes different of that of [19] - - for moderate and small chemical poten- 
tials. The former case might be applicable to the description of an epitaxial 
graphene subject to the gate voltage, while the later one is relevant for sus- 
pended graphene samples or for decoupled graphene layers [54]. 

Analyzing the conductivities similarly to the case of large /x, one can show 
that for small or moderate chemical potential the biggest contribution to the 
sum (42) in the vicinity of (43) comes from the no-term. However, in this 
case the extremum of the real part of G no /po lie at po ~ (M no+ i — M no ) ± T, 
where 

here we assumed that Y <C M no+1 — M no . For a small mass gap, A, this 
value is almost independent of the Landau level number, and thus the max- 
imum of rotation angle is inversely proportional to the difference between 
corresponding Landau levels 

UA-0)~T ^'t". (47) 

T y/n + 1 - Jn^ 




,«=170meV 



/J-lOmeV 



Figure 3: 9 and T as a function of frequency po[ me V] for different chemical 
potentials for B = 7T, V = lmeV, A = 0. 
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For A = also M = 0, so that two points, po — Mi + M and po = Mi — M 
coincide. Consequently, one has two identical contributions for n = 0, which 
is taken into account by 8o no in the formula above. 

The expression (47) gives a rough order of magnitude estimate of the 
Faraday rotation at the a 1 level. For describing the effect in very clean 
graphene samples (such as, e.g., reported in [54] where T ~ 3 • 10 _2 meV) 
one has to use complete expression (23) to correctly take into account the 
0(a/T) terms. 

For a suspended graphene we use (23) to plot on Fig. 3 both the Faraday 
rotation and transmission spectra against the frequency at B = 7T and 
r = lmeV for three different values of chemical potential. The maximum of 
rotation angle for such parameters has the order of magnitude predicted by 
(47), however, the minima in transmission spectra calculated via (23) are at 
least two times shallower than the ones predicted by (21). Still, both effects 
in this case are very large for an one-atom thick material, and presumably 
can be measured after minor modifications of the existing experimental set- 
ups. On the other hand, a very promising experiment was recently proposed 
in [20] . It was predicted there that with graphene placed in an optical cavity 
the Faraday rotation can reach huge values in the infrared region (up to 60 
degrees) and some modest values in the visible one. 

Finally, we can consider the rotation angle (21) as a function of chemical 
potential for a fixed frequency. One would naturally expect that for a small 
T its behavior should be similar to the one in the ideal case. Indeed, a 
numerical analysis shows that it is quantized in the same manner as the Hall 
conductivity as defined by (31). On the upper picture of Fig. 4 we plot the 
Faraday rotation fixing the frequency close to Mi. Thus for \x < Mi the 
effect gets enhanced according to (47). On the lower picture the frequency 
is far away from any special value, and thus the overall amplitude is much 
smaller. 



6 Conclusions 

In the present paper we recalculated the polarization operator of the Dirac 
quasi-particles in graphene in external magnetic field and established a clear 
connection between its components and optical properties of both suspended 
and epitaxial graphene. 

We showed, that in a number of different regimes the rotation of polar- 
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ization of light passing through monolayer graphene samples is giant being 
of the order of 0.1 rad, while the absorption can reach 40%. 

In particular we showed that the data of the recent experiment on Faraday 
rotation in the cyclotron resonance regime can be nicely fitted by predictions 
of the Dirac model for sufficiently high magnetic fields (starting with 2T). 
Moreover, we predict that for smaller chemical potentials (characteristic for 
suspended graphene) the effect is still very pronounced. We envisage the 
pole-like peaks in the polarization rotation angle and the transmission spec- 
tra at relatively high frequencies. We also reveal that in clean graphene 
samples the Faraday rotation as a function of chemical potential should be 
quantized in a similar fashion as the unconventional Integer Hall effect in 
graphene. 
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